tab br_lid_reg, nol
gen braslider=1 if  br_lid_reg==1
replace braslider=0 if braslider==.
gen faXUI_enc=gobop*UI_enc
gen ppnXUI_enc=ppn*UI_enc
gen ppcXUI_enc=ppc*UI_enc
set more off
gen mercosur=1 if cual_reg==3
replace mercosur=0 if mercosur==.

*modelo 0
logit tlc_china  mercosur inf_merc braslider, vce(r)
mfx
estimates store m0
esttab m0, aic

*modelo 01
logit tlc_china UI_enc  mercosur inf_merc braslider, vce(r)
mfx
estimates store m01
esttab m01, aic

*modelo 1
logit tlc_china UI_enc gobop mercosur inf_merc braslider, vce(r)
mfx
estimates store m1
esttab m1, aic

*modelo 2
logit tlc_china UI_enc ppn ppc reg mercosur inf_merc braslider, vce(r)
mfx
estimates store m2 
esttab m2, aic

*modelo 3

logit tlc_china UI_enc ppn ppc reg ppcXUI_enc ppnXUI_enc mercosur inf_merc braslider, vce(r)
mfx
estimates store m3

esttab m3, aic

*modelo 4
logit tlc_china UI_enc ppn ppc reg ppcXUI_enc ppnXUI_enc mercosur inf_merc braslider  rel_chi_hoy, vce(r)
mfx
estimates store m4
esttab m4, aic

*http://repec.org/bocode/e/estout/advanced.html#advanced001 
esttab m0 m1 m2 m3 m4,  pr2 aic se margin


*Figure 1
graph box inf_br inf_china inf_eeuu inf_eu , scheme(s2mono)

*Figure 2
prgen inf_merc, from(0) to(10) generate(prideo) rest(mean) ci  
graph twoway connected pridep1  pridep1lb pridep1ub pridex , ytitle("Pr. Y=1") xtitle("Mercosur's influence") scheme(sj) 

*Figure 3
logit tlc_china  i.ppn##c.UI_enc i.ppc##c.UI_enc reg mercosur inf_merc braslider, vce(r)
margins ppn ppc, at(UI_enc=(1(1)10))
marginsplot

